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Abstract 

We evaluate both the tau lepton energy loss produced by photonuclear inter- 
actions and the neutrino charged current cross section at ultra-high energies, 
relevant to neutrino bounds with Earth-skimming tau neutrinos, using differ- 
ent theoretical and phenomenological models for nucleon and nucleus struc- 
ture functions. The theoretical uncertainty is estimated by taking different 
extrapolations of the structure function F2 to very low values of x, in the low 
and moderate range for the tau lepton interaction and at high for the 
neutrino-nucleus inelastic cross section. It is at these extremely low values of x 
where nuclear shadowing and parton saturation effects are unknown and could 
be stronger than usually considered. For tau and neutrino energies E = 10^ 
GeV we find uncertainties of a factor 4 for the tau energy loss and of a factor 
2 for the charged current neutrino-nucleus cross section. 
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1 Introduction 



The detection of high energy neutrinos is one of the most important challenges in 
Astroparticle Physics. Conventional neutrino detectors exploit the long range of muons 
produced by muon neutrino charged current (CC) interactions [1]. With the discovery of 
neutrino flavor oscillations it has been realized that also tau neutrinos reach the Earth in 
spite of being heavily suppressed in all postulated production mechanisms. The possibility 
to search for tau neutrinos by looking for tau leptons that exit the Earth, Earth- skimming 
neutrinos, has been shown to be particularly advantageous to detect neutrinos of energies 
in the EeV range [2, 3]. The short lifetime of the tau lepton originated in the neutrino 
charged current interaction allows the tau to decay in flight while still close to the Earth 
surface producing an outcoming air shower in principle detectable by both fluorescence 
telescopes and air shower arrays [4]. This same channel yields negligible contributions 
for other neutrino flavors. The sensitivity to tau neutrinos through the Earth-skimming 
channel directly depends both on the neutrino charged current cross section and on the tau 
range (the energy loss) which determine the amount of matter with which the neutrino 
has to interact to produce an emerging tau [5, 6]. While the energy loss for muons is 
shared by roughly equivalent contributions from pair production, bremsstrahlung and 
photonuclear interactions, for tau leptons of energies above E = 10^ GeV, photonuclear 
interactions (i.e. lepton-nucleus inelastic interactions dominated by small values of Q^) 
are responsible for the largest and the most uncertain contribution [7, 8, 9] . 

Both the neutrino cross section and the tau photonuclear energy loss are calculated 
from theory using structure functions which carry the information of the nucleon and nu- 
cleus structure. In order to study the uncertainties in the calculation of Earth-skimming 
neutrinos the same structure functions should be consistently used for both processes due 
to their strong correlation in the resulting tau flux. Unfortunately this is not possible 
since the kinematical (minus the squared momentum transfer) and Bjorken-x ranges 
that contribute to these processes are quite different, specially at EeV energies, and the 
available parameterizations are not entirely adequate to describe both ranges simultane- 
ously. 

The scale that contributes to the tau energy loss, dominated by photon exchange. 
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is low and moderate at very low where perturbative and non perturbative QCD 
effects are mixed. The CC neutrino cross section is produced by IV-boson exchange that 
sets the relevant scale of Q'^ to values up to at low a region where perturbative 
QCD is expected to work. In both cases the relevant x range hes well outside the regions 
where structure functions are measured, so one has to rely on extrapolations which contain 
significant uncertainties. 

The charged current neutrino cross section is usually calculated using parton distribu- 
tion functions which are evolved according to perturbative QCD predictions. A number of 
alternative parameterizations exist, some of which allow extrapolation of the uncertainties 
in the fitted parameters as a mean to explore some of the uncertainties associated to the 
calculation. In the case of photonuclear processes existing predictions at high energy arise 
basically from two independent approaches, the Generalized Vector Dominance (GVD) 
model and Regge-hke models. 

In this article we study the tau energy loss (see also Ref. [10]) and the neutrino- 
nucleus cross section. Both quantities have direct implications for high energy neutrino 
detection, in particular for Earth-skimming tau neutrinos. Due to the large uncertainties 
in the existing models, the fact that none of them covers simultaneously the kinematical 
region relevant for both quantities, and the need of consistency in both calculations, we 
use and extend available models with the aim of estimating the theoretical uncertainty 
by considering extreme results. In this way, in the frame of the most relevant models, 
we cover the range of possible scenarios for the extrapolation of structure functions to 
the relevant x and range. Two important effects to be taken into account in this 
extrapolation of the structure functions are nuclear shadowing corrections and saturation 
due to partonic screening. Nuclear corrections [11] are deviations from the naive picture in 
which the nucleus is treated as an incoherent sum of nuclcons. Saturation [12] accounts for 
the fact that the structure functions cannot rise indefinitely as x goes to zero. Saturation 
effects may be included in nuclear corrections but are also present in the nucleon structure 
functions (although for smaller values of x and/or Q^). In addition to existing calculations, 
a new computation of the tau energy loss and the neutrino-nucleon cross section based 
on saturation physics [13] is also presented in this work. 

The result of the present analysis is an uncertainty band for both the tau-lepton 
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energy loss by photonuclear interactions and the CC neutrino-nucleus cross section. Un- 
derstanding and minimizing the uncertainties in these two calculations must be considered 
an important priority for high energy neutrino astrophysics. 



2 The photonuclear tau energy loss 

The average energy loss per unit depth, X, of taus is conveniently represented by: 

-(^^^ = a{E) + b{E)E, (1) 

where a{E) is due to ionization and b{E) is the sum of fractional losses due to e+e~ 
pair production, bremsstrahlung, and photonuclear interactions. The parameter a{E) is 
nearly constant and the term h{E)E dominates the energy loss above a critical energy 
that for tau leptons is of a few TeV. The electromagnetic contribution to the energy 
loss, mainly due to pair production and bremsstrahlung, is well under control, while the 
photonuclear interaction which dominates for tau energies exceeding E = 10^ GeV is 
affected by relatively large uncertainties. 

The contribution to h{E) from photonuclear interactions is obtained by integration of 
the lepton-nucleus differential cross section, da^^/dy: 

where Na is Avogadro's number, A the mass number, and y the fraction of energy lost 
by the lepton in the interaction. 

For the lepton-nucleus differential cross section we consider the general expression for 
virtual photon exchange in terms of structure functions: 



dQ^dy y 



(3) 



AE'^ \ 2(1 + R^) 

where E is the lepton energy in the lab frame, mi the lepton mass, and a the fine 
structure constant. is the structure function F2 for a nuclear target A which is 
found to be different from the mere superposition of A free nucleon structure functions 
F2 [11]. R"^ is the ratio of the longitudinal to transverse structure functions which gives 
a small contribution to the cross section [7] and is neglected for clarity of the discussion 
below. The variables x, y and are related by kinematics through = 2MExy, and 



both F2 and R arc functions of x and . The contribution to the tau energy loss from 
neutral current and 7-Z interference interactions was estimated to be small [14] and is 
also neglected. 

The hmits in the double integral of Eq. (2) are well estabhshed: 

2 2 

Qmin = ' ^"^"^ " 2mpEy - 2m„mp - ml , (4) 

_ 2m^iii^, + ml - -I _!!}! /'c:^, 

Vmin ^ ' Umax ^ ' 

where and m,r are the proton and pion mass, respectively. 

The predictions of the photonuclear interaction cross section in the GVD Model [15] 
(BB) and in its extension to higher energies by including a perturbative component based 
on the color dipole model [16] (BS), have been widely used to explore muon and tau lepton 
propagation in matter (see for instance [17, 8, 9] and references therein). 

The calculations in which the F2 structure function is given by a phenomenological 
parameterization of data based on Regge Theory appear in Refs. [7] (DRSS), [14] (BM), 
[18] (KLS), and [19] (PT). For the proton structure function, F2, DRSS (see also Ref. 
[20]) uses the ALLM model [21], while BM and KLS both consider the CKMT model [22] 
at low matched at high to perturbative QCD predictions based on different param- 
eterizations of parton distribution functions, and PT uses the proton structure function 
of Ref. [23]. The structure function is shown in Figs. 1 and 2, together with the HERA 
data at the lowest measured x values at different Q^. 

In DRSS, BM, and KLS calculations the nuclear structure function is related to the 
proton structure function through — /"^AFf . At high energy only the low x behavior 
of the nuclear correction factor is relevant to the calculation of b{E), as we will show 
below (see Fig. 6). In the DRSS calculation the low x behavior of freezes at the value 
yA ^ ^-0.1 fQj. ^ ^ 0.0014 {r^ 0.73 for standard rock, ^ = 22), while in the BM (and KLS) 
calculations reaches a maximal asymptotic regime = A~^/^ (~ 0.36 for A — 22) at 
much lower x (see Fig. 3). Both DRSS and BM nuclear corrections are Q^-independent. 

In addition to the existing calculations we present a new computation of the pho- 
tonuclear tau energy loss using the results of Ref. [13] (ASW) which arc based on the 
geometric scaling property [26] that all data on a'^*^ and on cr'"'*^ lie on a single universal 
curve in terms of the scaling variable r = Q'^ iQlat whose form is inspired in saturation 



physics (the detailed expressions leading to the ASW F2 structure function are given in 
the Appendix). The ASW F2 structure function for the proton case is plotted in Figs. 1 
and 2 (for x < 0.01 where this parameterization is expected to be vahd). The ASW struc- 
ture function F2 contains mild nuclear corrections at low x when compared with DRSS 
and BM nuclear corrections (see Fig. 3). Nuclear corrections in ASW depend on Q^. 

The photonuclear contributions to h{E) computed (for standard rock A = 22 through- 
out all this paper) with ALLM and with CKMT structure functions, and the same nuclear 
corrections [7], give very close results (see Fig. 4). Although ALLM and CKMT parame- 
terizations share a common theoretical base, with a reggeon and a pomeron component, 
and they are fitted to the same data sets, ALLM systematically hes above CKMT at low 
X (see Fig. 2), which accounts for the difference in h{E) observed in Fig. 4. 

The lowest values of h{E) at high energies is obtained with the ASW structure func- 
tions. Though the ASW structure function F2 contains mild nuclear corrections at low 
X, saturation effects at the nucleon level are rather strong and limit the rise of h{E) 
with energy as observed in Fig. 4. For energies below E = 10^ GeV the result from the 
ASW structure function is higher than those from ALLM or CKMT (see Fig. 4). This 
is because at low the ASW structure function is significantly higher for the region 
10"^ < X < 10~^ (see Fig. 2) which is the relevant range for energies below E — 10^ GeV, 
as it can be deduced from Fig. 6. Thus the saturation-based ASW prediction lowers the 
energy loss rate h{E) with respect to the already existing predictions by a factor 2 at 
E = 10^ GeV, and by a factor even larger at higher energies. 

The BB/BS calculation gives the largest of the predicted energy loss rates up to 
energies of the order E — 10^ GeV. Above this scale the PT result exceeds all other 
existing predictions by at least a factor 2 already at £^ = 10^ GeV (i.e. a factor 4 with 
respect the ASW prediction, see Fig. 4). Thus the PT prediction can be considered as an 
estimate of the upper limit of the tau energy loss at UHE. 

Much of the uncertainty in the tau energy loss is actually due to nuclear effects. The 
choice of nuclear corrections from Ref. [7], Ref. [14], or from Ref. [13] (see Fig. 3), 
translates into differences in the calculated value of b{E) (using the ALLM structure 
function) by a factor rising from 1.5 to 2.5 as the tau energy increases in the range 
E — 10^-10^ GeV (see Fig. 5). This energy range corresponds to the region of very low 
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X where differences in the nuclear correction factor are large. In order to quantify how 
much different regions of x and contribute to b{E), the dependence of b{E) on the 
maximum value of x and on the maximum value of considered in the integration is 
shown in Figs. 6 and 7. 

The differential cross section da'^'^/ dy is also a relevant quantity for high energy neu- 
trino detection as it enters the event rate convolutions together with the neutrino flux 
and the experimental acceptances. Indeed it has been shown that stochastic effects of 
the tau energy loss distribution have significant relevance in the prediction of emerging 
tau rates [20]. The energy loss spectrum yda'^^/dy obtained using both ALLM and ASW 
structure functions are compared in Fig. 8. Clearly, the energy loss spectrum calculated 
with ALLM is significantly harder than the one calculated with ASW. The contributions 
of moderate (Q^ > 1) and low Q"^ {Q^ < 1) (in a rough way corresponding respectively 
to hard and soft interactions) are shown separately in Fig. 9 for the ALLM structure 
function. 

3 The charged current neutrino cross section 

The absolute value of the cross section naturally has a direct impact on the sensitivity 
of experiments because the event rate is directly proportional to it, but it also enters with 
opposite effect in the attenuation of the neutrino beam as a function of matter depth 
traversed, having much impact on the angular and energy distribution of the events. These 
two effects combine in the case of Earth-skimming tau neutrino interactions to play an 
important role for the rate calculation. In addition to the tau lepton photonuclear cross 
section we also study how the uncertainties in the F2 structure function at low x affect 
the CC neutrino deep inelastic cross section. Since in the more realistic expectations [27] 
the nuclear corrections to the CC neutrino-nucleon cross section decrease at low x with 
increasing Q^, becoming small at high Q'^ [28], we will neglect them in our calculations. 

The CC DIS neutrino-nucleon cross section is expressed in terms of the structure 
function F2 as follows: 



where E is the neutrino energy and y the fraction of energy lost by the neutrino in the 
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interaction. In this expression and xF^ contributions are neglected since Fl tends to 
zero as Q"^ rises and xF^ deals basically with the valence partons which hardly contribute 
at the low x values relevant for the cross section. 

In order to consistenly use the structure functions from charged lepton interactions 
(as ALLM, CKMT, and ASW models) in neutrino interactions we must relate the electro- 
magnetic and weak structure functions. The F2 structure function for neutrino interaction 
is related to the F2 structure function for charged lepton interactions by the ratio of the 
weak and electromagnetic couplings through = 18/5 (assuming a symmetric 
sea) , although the kinematical regions of the two processes are different (low and moder- 
ate Q"^ ~ 0.01-10 GeV^ in the photonuclear case and high ~ in the high energy 
CC interaction). Concerning the x range the main contribution comes from low x in both 
cases, though x values are lower in the photonuclear case than in the CC interaction. For 
the calculation of the neutrino-nucleon cross section at high energies, we then use the 
structure function F2 for charged lepton interaction valid up to very low x and high Q^, 
instead of following the standard approach based on parton densities. 

The neutrino-nucleon cross sections from ALLM and CKMT structure functions are 
presented in Fig. 10. They are clearly below predictions from modern parton densities 
[29], since the ALLM parameterization is not consistent with high experimental points 
(see Fig. 1) and CKMT is not evolved to high Q^, so we have not used them to discuss 
the theoretical uncertainties in the estimation of the CC neutrino-nucleon cross section. 
Instead we have taken the parameterization of F2 a la, BCDMS obtained by the SMC 
Collaboration [30], which correctly represents the existing experimental data at high 
(see Fig. 1) and provides a smooth connection at neutrino energies around E — 10^ GeV 
with the parton density prediction of the CC neutrino-nucleon cross section [29] (see 
Fig. 10). 

We have performed three different extrapolations at low x of the F2 parameterization 
given in Ref. [30], one following the ASW structure function, a second one from the 
phenomenological parameterization fitting low x HERA data [31], and the third one which 
corresponds to the double logarithmic approximation (DLA) in QCD from Ref. [32] 
(KOPA) (ASW and KOPA structure functions are valid at low x, x < 0.01, i.e. at high 
energy) . 
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In Fig. 11 the effect of taking the three different parameterizations of the structure 
function F2 at low x on the neutrino-nucleon cross section is shown. We see that in 
comparison with the prediction obtained with evolved QCD parton densities of Ref . [29] , 
both KOPA (which corresponds to the DLA of perturbative QCD) and ASW (which 
includes strong saturation effects) estimations are below at high energies. 

On the other hand the extrapolation of the HERA based parameterization with the 
exponent A — 0.0481 ln((5^/0. 292^) (F2 ~ a^~^), produces an extremely fast increase 
of the cross section with energy (see the upper curve in Fig. 11), since this exponent 
rises to values above A ~ 0.5 when Q'^ becomes large. This raw extrapolation is in 
contradition with perturbative calculations and we do not consider it for uncertainty 
estimates as it is not physically motivated. Nevertheless it is considered here to explicitely 
show its discrepancy with pQCD. For the more realistic scenarios, when the rise of the 
exponent freezes to smaller values A < 0.4, our prediction supports the result obtained 
in the detailed analysis of Ref. [29]. The curves are shown in Fig. 11 (from up to down 
corresponding to A frozen to A = 0.50, 0.40, and 0.38 respectively). When considering 
only physically motivated extrapolations, the theoretical uncertainty at E = 10^ GeV is 
a factor 2. 

4 Conclusions 

We estimate the uncertainties coming from the extrapolations of the existing models 
for proton and nucleus structure functions for tau energy loss and for CC neutrino-nucleon 
cross section. Both calculations must be done consistently within the same model as their 
effect on the tau flux produced by Earth-skimming neutrinos is correlated. The theoretical 
uncertainty in the tau energy loss is greater than that of the neutrino-nucleon CC cross 
section because the region contributing to the tau energy loss cross section is lower and 
so are the relevant values of Bjorken-x. In addition the structure functions conventionally 
used for the calculation of the tau energy loss are not suitable to be used in the high-Q^ 
range which is relevant for the CC neutrino-nucleon interaction. As a result systematic 
effects arising in the calculation of a tau neutrino bound from Earth-skimming events due 
to uncertainties in the structure functions turn out to be difficult to evaluate. Several 
extreme models allowed by extrapolation of structure functions have been explored in 
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order to estimate ranges for these quantities. 

Below energies in the E — 10^ GeV range the highest prediction for the photonuclear 
contribution to tau energy loss, b{E), is provided by the BB/BS calculation. Above this 
energy range the PT result exceeds all other considered predictions while the lowest cal- 
culation is obtained using the ASW structure functions. The difference between the two 
extreme predictions reaches a factor 4 at ii^ = 10^ GeV and increases as the energy rises. 
The BB/BS, ALLM, and CKMT calculations agree within a 30 % and go approximately 
parallel for all energies, which is an indication of a systematic normalization difference of 
the structure functions in each model. The application of much stronger nuclear shadow- 
ing (than usually considered) at low x can lower the prediction of b{E) with respect to 
the already existing calculations by a factor up to 2 at = 10^ GeV. 

In the case of the CC neutrino-nucleon cross section the importance of nuclear effects 
at high energies is expected to be small [28] . We have also considered saturation effects in 
the CC neutrino-nucleon cross section by using the structure function ASW. At £■ = 10^° 
GeV, the CC neutrino-nucleon cross section calculated with the ASW structure function 
is found to be half of the pQCD calculation with parton densities, in rough agreement 
with the evaluation of saturation effects reported in Ref. [33], and also in Ref. [34] where 
different parameterizations of the dipole cross section containing saturation are employed. 
The calculation of the neutrino cross section with the ASW structure function has also 
been performed in Ref. [35]. Though the quantitative agreement of our result with this 
calculation is reasonably good, some discrepancy appears due to the fact that to extend 
the validity of ASW to the region x > 0.01, we have connected the ASW structure function 
to the parameterization of HERA data d la BCDMS from Ref. [30] . 

The effect of a rapid rise of the F2 structure function at low x in the CC neutrino- 
nucleon cross section has also been studied using the x-slope A(Q^) of the F2 HERA data 
{F2 ~ a;"'^) for all values. We have found that the cross section rises with energy very 
rapidly. At = 10^° GeV it can become a factor 4 above the pQCD calculation with par- 
ton densities. On the other hand the logarithmic rise of the structure function F2 at small 
X predicted by the DLA-pQCD results in a slower increase of the CC neutrino-nucleon 
cross section with energy. At E = 10^^ GeV the DLA estimation is a 20 % below the 
pQCD calculation with parton densities. When considering only realistic extrapolations. 



9 



the theoretical uncertainty aX E = 10^ GeV is a factor 2. 

The obtained uncertainty for the tau energy loss is to be implemented, together with 
the corresponding one for the CC neutrino-nucleon cross section, both in analytical and 
Monte Carlo calculations of the rates of taus emerging from Earth-skimming tau neutrinos, 
which is currently being used to calculate high energy neutrino bounds. This task is 
beyond the scope of this paper. 

Acknowledgements 

We thank O. Blanch Bigas, M.V.T. Machado, D. Pertermann, Yu.M. Shabelski, and 
D.A. Timashkov for useful comments on this work. N.A. acknowledges financial support 
by Ministerio de Educacion y Ciencia (MEC) of Spain under a Ramon y Cajal contract. 
This work has been supported by MEC under grants FPA2005-01963 and FPA2004-01198, 
by Xunta de Gahcia under grant 2005 PXIC20604PN and Conselleria de Educacion, and 
by FEDER Funds. 



10 



Appendix: The ASW F2 structure function 



The form of the single universal curve where all data on a"^*^ and on a^*^ lie as function 
of the scaling variable r = /Qlat motivated by saturation and given by [13, 36]: 

<7^>(x, Q^) = $(r) = ao [^e + T (0, i) + In^] , (7) 

with ^jE the Euler constant, V (0,,^) the incomplete F function, and ^ = a/r^, with a = 
1.868 and b = 0.746 extracted from a fit to lepton- proton data. The saturation scale Q^^^ 
is parameterized as Ql^^{GeV^) = (x/xq)-^ [37], where xq = 3.04 ■ 10"^, A = 0.288, and 
X — X {Q^ + Am'f)/Q^ with nif — 0.14 GeV. The normahzation is fixed by CTq = 40.56 //b. 
The extension to the nuclear case is done through 

= (8) 

and 

QU. - QU. (^) - - [^] . (s) 

where the nuclear radius is given by the usual parameterization Ra — (1.12^4^/^ — 
0.86A-V3) 5 ^ 0.79 ± 0.02 and nR^ = 1.55 ± 0.02 fm^ are extracted from a fit to 

Icpton-nucleus data. The nuclear structure function is F^^XjQ"^) = Q"^ a"'* / {A-n"^ a) . 
The ASW structure function for the proton case is recovered by taking ^4 = 1 in the 
expressions above (see Figs. 1 and 2). 

The functional shape of (7) is motivated by considerations in saturation physics [13, 
36] . From a pragmatic point of view, it provides a very good description of existing lepton- 
proton and lepton-nucleus data in the region 0.01 < r, ta < 100 and x < 0.01 which for 

= 0.01,0.1, 1, and 10 GeV^ corresponds to a low x limit of ~ 10-^ 10"^ 10"^°, and 
10~^^, respectively. For r ^ F2 behaves like a single logarithm, so F2 oc Inl/x for 
x — > and F2/A oc In 74/^4^/^ for ^4 — > 00. Thus this model results in very large screening 
corrections for asymptotic values of x and A. 
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Figure 3: The nuclear correction factor as a function of x. 
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Figure 4: The photonuclear energy loss rate, b{E), computed in different models. 
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Figure 5: The effect of nuclear corrections on the photonuclear energy loss rate, b{E). 
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Figure 8: Spectrum of the tau energy loss by photonuclear interactions for tau energies 
E = 10^ 10^, and lO^^ GeV. 
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Figure 9: Spectrum of the tau energy loss by photonuclear interactions for a tau of energy 
E = 10^^ GeV. The contributions from values above and below = 1 GeV^ are shown 
separately. 
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Figure 10: The neutrino-nucleon CC cross section as a function of the neutrino energy, E, 
from: ALLM, CKMT, SMC d la BCDMS structure functions, and the parton density 
calculation by Anchordoqui et al. 
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Figure 11: The neutrino-nucleon CC cross section as a function of the neutrino energy, E, 
from: ASW (extrapolation based on saturation physics), KOPA (extrapolation based on 
DLA QCD) , HERA (phenomenological parameterization of HERA data) , and Anchordo- 
qui et al. (parton density calculation). 
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